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We review the methods and results obtained in an analysis of the experimental heavy ion 
collision research program at nuclear beam energy of 160-200^4 GeV. We study strange, 
and more generally, hadronic particle production experimental data. We discuss present 
expectations concerning how these observables will perform at other collision energies. 
We also present the dynamical theory of strangeness production and apply it to show 
that it agrees with available experimental results. We describe strange hadron production 
from the baryon-poor quark-gluon phase formed at much higher reaction energies, where 
the abundance of strange baryons and antibaryons exceeds that of nonstrange baryons 
and antibaryons. 
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1. Introduction and Overview 
1.1. Introduction 

The study of highly excited and dense hadronic matter by means of ultra- relativistic 
nuclear collisions is a relatively novel area of research at the border between nuclear 
and particle physics. As such it, is in a rapid experimental and theoretical evolu- 
tion. The primary goal of research in this field is the creation and investigation 
of elementary (particle) matter under extreme density and temperature conditions. 
The existence of a novel non-nuclear high temperature phase of elementary matter 
is an unavoidable consequence of the current knowledge about the strong nuclear 
interaction, rooted in the theory of strong interactions, the quantum field theory of 
quarks and gluons called quantum chromodynamics (QCD).EJu 

Discovery and study of quark-gluon plasma (QGP), a 'deconfined' state consist- 
ing of mobile, color-charged quarks and gluons, is the objective of the relativistic 
heavy ion experimental research program underway at the Relativistic Heavy Ion 
Collider (RHIC) at Brookhaven National Laboratory (BNL), New York, and at the 
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Super Proton Synchrotron (SPS) accelerator at the European Organization for Nu- 
clear Research (CERN), Geneva.0 In a recent half day long workshop, and in an 
accompanying press release in February 2000 the CERN laboratory has formally 
announced that it views the collective evidence available today in their seven rela- 
tivistic nuclear collision experiments as being conclusive: ". . . A common assessment 
of the collected data leads us to conclude that we now have compelling evidence that 
a new state of matter has indeed been created, at energy densities which had never 
been reached over appreciable volumes in laboratory experiments before and which 
exceed by more than a factor 20 that of normal nuclear matter. The new state of 
matter found in heavy ion collisions at the SPS features many of the characteristics 
of the theoretically predicted quark-gluon plasma. ..." .0 

This research program has been developed over the past two decades in order 
to study the properties of elementary matter at conditions similar to those seen in 
the very early universe 30^s after the big-bang, before the temperature decreases to 
about T = 150MeV~ 1.7- 10 12 K. It has been many times argued that it is possible to 
achieve a laboratory recreation of this condition in a small-bang relativistic nuclear 
collision. The question is at what collision energy the transition to a color deconfmcd 
QGP phase first occurs. Early suggestion has been that this could occur at an 
intrinsic available energy per participating nucleon as low as 4 and 8 times the 
nucleon mass, corresponding^ to the range 30 and 120 A GeV per nucleon beam energy 
in fixed target experiments.u 

The conditions in the early universe and those created in nuclear collision ex- 
periments differ somewhat: whereas the primordial quark-gluon plasma survived 
for about 30 /is in the big bang, the comparable conditions in nuclear collisions are 
not expected to last for more than 10 -22 s due to the rapid explosion of the hot 
matter "fireball" . Moreover, in the matter created in heavy ion collisions quarks 
are expected to outnumber antiquarks noticeably due to the baryon content of the 
colliding nuclei, whereas the net relative excess of the quarks over antiquarks in the 
universe was less than 10 -9 . 

Numerical simulations of QCD suggest that the nature of the transformation 
between the hadronic and quark-gluon phases can change drastically as the values 
of the parameters of the theory are varied. Q'Q'H Recent analytical studies of the phase 
properties of QCD have supported the conclusion that the dependence on the net 
baryon density (baryochemical potential) is especially interesting.BE3 Perhaps the 
most fundamentally important observable in this context is the latent heat associ- 
ated with the breaking of color bonds among quarks, leading to the deconfinement 
of quarks. An experimental determination of this quantity and its dependence on 
beam energy would be of great scientific interest. 

Strange particle signatures for the formation and evolution of the deconfined 
quark-gluon phase of elementary matter form a significant cornerstone of experi- 
mental QGP discovery. This subject has been developed quite intensely for the past 
20 years. EiESll^'ESEjtJ The enabling difference in physics between confined and de- 
confined matter concerning strange particle signatures is rather simple: 
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In the QGP phase the particle density is high enough and the strange flavor 



ce of stran- 
whilc in the 



production energy threshold low enough to assure that a high ahunda,n < 
geness can actually be produced on the time scale available^llilij£3ll3 \ 
confined matter phase this has been shown not to be the casejiM2l as long as one 
wants to remain consistent with other experimental results. 

• Population at, and even in excess, of chemical equilibrium of hadron phase space 
occupancies occurs only when the entropy rich QGP phase disintegrates rapidly and 
explosively into hadrons ME3ii3'E3 

There are several important and when viewed together, uniquely QGP charac- 
teristic predictions regarding strangeness, expected to occur should deconfincmcnt 
set in. Specifically, the three pillars on which the QGP hypothesis stands when seen 
by means of strangeness flavor observables are: 

1) matter-antimatter symmetry as seen for directly emitted strange baryon and 
antibaryon particles in the m^-spectral shape and strange quark fugacity; 

2) (multi)strange baryon and antibaryon enhancement increasing with strangen- 
ess content; 

3) enhancement of the (specific) strangeness flavor yield per reaction partici- 
pant (baryon), by a factor 1.5-3 at SPS conditions, the value depending on what 
is used as baseline, and if one looks alone at the central rapidity region, where this 
effect is strongest, or considers the global strangeness yield, including the kinematic 
domains of projectile and target fragmentation. 

All three predictions have recently been confirmed at the current SPS energy 
range 158 GeV per nucleon for Lead (Pb), and some also for 200 GeV per nucleon 
Sulphur (S) induced reactions. 
Examples and stepping stones are in particular: 

1) The WA97 collaboration reported a detailed study of transverse mass strange 
baryon and antibaryon spectra which_show a highly unusual symmetry between 
strange baryon and antibaryon sector.Eil 

2) A detailed analysis of Pb-Pb results by the WA97 collaboration has demon- 
strated, comparing p-p, p-A with A-A results, a strong enhancement in the per- 
tinent (rrmlti ) s tr a n g e baryon and antibaryon yields, increasing with strangen- 
ess content .ESfpHil'Ea The results of the NA49 collaboration are consistent with 
these findings.ES The WA85 collaboration also finds an enhancement of multi- 
strange baxyons and antibaryons, increasing with strangeness content in S-induced 
reactions. Ell 

3a) Strangeness enhancement at, m id-rapidity has been observed imS-induced re- 
actions by the experiments NA35,N WA85 and WA94|!3 and NA44E9 In thelarffer 
Pb-Pbreaction system strangeness enhancements are reported by WA97,c3t3H 
NA49,E3 and NA44.E3 Results of the experiment NA52 suggest further that the 
onset of strangeness enhancement occurs rather suddenly as the centrality of the 
collisions and thus the size of participating matter rises above baryon number B — 
40-50 .0 

3b) Global Strangeness enhancement has been observed both in S-induced and 
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in Pb-Pb reactions by the NA35,EI and NA49J 3 ! experiments. 

In this article, we also rely in many aspects of this discussion indirectly on other 
experimental results of NA35 and NA49 collaborations, which offer a globalview 
on particle production pattern considering the large kinematic acceptance P 5 " 3 ^ We 
will not discuss or use here other experimental discoveries which have contributed 
to the CERN announcement, which are not related to strangeness, such as J/*f> 
suppression, dilepton and direct photon production. Readers interested in these 
topics should consult the list of experimental resultso 

In view of many often intricate but, when analyzed, convincing experimental 
findings about strange particle production, the purpose of this article is to present 
a comprehensive and selfconsistent view on the understanding of the evidence com- 
prised in strange hadron production for the formation of quark-gluon plasma at 
CERN, and to discuss resulting expectations how this observable will perform at 
RHIC. We address in this review: 

i) in section 2 the status of a analysis of the experimental data obtained at SPSJ23 

ii) in section 3 the implications of these results for the understanding of the dense 
phase formed in these reactions, 

iii) in section 4 an adaptation of the dynamical theory of strangeness production 

in QGP to RHIC conditions, 

iv) in section 5 an application of these results to obtain predictions for hyperon 
yields from QGP at RHIC, 

v) in section 6 highlights of the results presented here and we draw our conclusions. 

1.2. Overview 

We introduce in section 2, the Femii-2000 model,0@ a straightforward elabora- 
tion of the original Fermi proposal,Ea that final state strongly interacting particles 
are produced with a probability commensurate to the size of the accessible phase 
space. In this approach, the hadron phase space is characterized to the required 
degree of accuracy by six parameters which have a clear physical meaning and can 
be in computed, and/or rather easily understood qualitatively, when their values 
have been determined by an analysis of the experimental data. The physical picture 
underlying the use of the statistical Fermi model in the 21st century is the sudden, 
explosive disintegration of a high temperature hadronic matter fireball, apparently 
consisting of deconfined quark-gluon matter. Since its proposal 50 years ago, the 
Fermi model has been subject to considerable scrutiny and adaptation, withHage- 
dorn's 'boiling' hadronic matter being the most important stepping stone.Ej The 
following were the relevant recent steps in the development of the statistical par- 
ticle production description required to analyze the strange particle experimental 
results: 

1. Considering the interest and considerable theoretical effort vested in under- 
standing strange quarknroduction mechanisms and the study of chemical 
equilibration processes,t3'Ej it was a natural refinement to introduce an ex- 
pression of chemical non-equilibrium in the number of strange quark pairs, 
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7 S 7^ into the analysis of strange hadrons. 

2. In the analysis of the entropy content in S-W/Pb interactionsjil we found 
entropy excess related to excess of meson abundance. This prompted us to ex= 
plore possible nonequilibrium yield of mesons compared to that of baryons.cd 

3. Since the dense hadron fireball is subject to explosive disintegration, final 
state hadrons emerge from rapidly outward moving volume cells. In order 
to describe quantitatively spectra of hadronic particles, and their yields in 
restricted domains of phase space, such 'collective' matter flow motion needs 
to be modeledSffi 

4. Chemical (particle abundance changing) processes occurring at the time of 
hadronization do not generally lead to an equilibrium chemical yield of light 
quarks, and the chemical non-equilibrium is more pronounced if hadronization 
is a sudden process on the time scale of chemical quark equilibration. While 
this effect has been well accepted for strange quarks, as these need to be 
produced in microscopic processes, the need to introduce the light quark pair 
abundance parameter, 7 g ^ 1, was recognized rather late,0cl and is not yet 
widely accepted. 

A data analysis we perform, allowing for all these effects, does not simply yield a 
set of 'best' parameters; rather: 

i) it offers a complete characterization of the phase space of hadrons and its 
occupancy, allowing one to extrapolate reliably the particle yields into kinematic 
domains not accessible at present; 

ii) it allows one to study and understand the magnitude of model parameters, so 
that we can safely to extrapolate their values to other reaction conditions, e.g., from 
SPS to RHIC as will be done here; 

iii) it allows one to evaluate the physical properties of the phase space characterized 
by these parameters, which provides very precise information about the physical 
properties of the particle source. This in turn leads to the understanding of the 
nature of the dense fireball created in the heavy ion collision. 

In order to pursue these aims, we need to reach considerable precision in the 
description of experimental results. 

Our approach and objectives elaborate significantly on the now commonly ac- 
cepted observation that all hadronic particles produced in strong interaction pro- 
cesses qualitatively satisfy statistical model predictions, as d iscovered and discussed 
in great detail by Rolf Hagedorn more than 30 years agoria The fajcLihat the statis- 
tical model indeed 'works' does not cease to amaze and impress. c£Ie3 At times this 
even provokes the hypothesis that 'thermal' abundances of hadrons could arise in 
some mysterious and unknown way,c3 and thus one could proceed to predict 'ther- 
mal' yields of hadrons, apparently believing that the statistical 'thermal' model 
substitutes for conventional particle production mechanisms. The recent proposal 



6 Jean Letessier and Johann Rafelski. 



by Bialas to consider the fluctuation of stringtension is indeed suggesting how such 
an explanation could arise in p-p reactions J13 

However, for reactions of large nuclei, the statistical description implies and 
exploits the result of repeated occurrence of microscopic collision reactions, and 
the associated approach to an equilibrium distribution shape, and independently, 
also approach to chemical particle abundance equilibrium. This is independent of 
the above described possibility that in elementary reaction systems equilibration is 
possibly a consequence of microscopic properties of strong interactions. Our dis- 
cussion of A-A reactions thus aims at an improvement of the statistical description 
beyond the 'thermal' model, such that we can deal with standard deviation errors 
comparing theory and experiment, as is common in the field of particle and nuclear 
physics. 

The way we set up the Fermi-2000 model represents the microscopic processes 
that are occurring, and of course there are limits to this description. For exam- 
ple, primary high energy initial interactions can produce heavy quark flavor which 
could not arise from the generally softer interactions occurring in the kinetically 
equilibrated system (an example is expected production of charm at RHICH'til). 
More generally, an acceptable failure of a statistical description is the one which 
under-predicts the yield of rarely produced particles. For this reason, it is necessary 
to scrutinize the validity of the statistical description, at least for the most rarely 
produced hadrons. 

Pertinent results of a complete analysis of the Pb-Pb system are presented in 
section 3. We describe abundances and spectra of hadronic particles observed by 
both the wide acceptance NA49-experiment and the central rapidity (multi)strange 
(anti)baryon WA97-experiment. Our method of analysis shows that both these 
families of results obtained with widely different methods are consistent and it al- 
lows us to reach the required precision in their description. This objective could 
be reached only after we have introduced light quark chemical non-equilibrium and 
allowed that the all strange fl(sss) and fl(sss) hadrons are enhanced beyond their 
statistical phase space yield, a point we will discuss in greater detail in subsec- 
tion 3.4 below. These developments occurred after the last comprehensive review 
of the subject appeared,E3 and after the extensive SPS-Pb-beam experimental re- 
sults became available. The introduction of light quark chemical nonequilibrium 
has had a significant impact on the determination of the physical properties of the 
hadron emitting source, as it allows for a considerable reduction of the chemical 
freeze-out temperature: specifically we have determined that Tf = 143 ± 5 MeV, 
we have also included here an estimate of the systematic error, for the temperature 
at which practically all strange hadrons are formed. This result is nearly 30 MeV 
below values th a t one might infer otherwise in a qualitative study of the statistical 
hadron yields .EJO'E^I 

This relatively low freeze-out temperature is consistent with the result that the 
chemical freeze-out parameters determine correctly the shape of hadron mj_-spectra, 
which suggests that after the deconfined QGP source has dissociated the resulting 
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hadrons are practically free-streaming and thus that thermal and chemical freeze- 
out do not di ffer m uch if at all. This particle production scheme is called sudden 
hadronizationr3'E3 In terms of a microscopic model it occurs when hadronic particles 
are produced either in: 

a) an evaporation process from a hot expanding surface, or 

b) a sudden global hadronization process of exploding, possibly super-cooled 
deconfmed matter. 

Experimental evidence supporting the picture of sudden QGP hadronization 
is most directly derived from the baryon-antibaryon transverse mass m^-spcctral 
symmetry. Another piece of evidence for a sudden hadronization is the chemical 
overabundance of light quark pairs in hadronization and the associated maximiza- 
tion of entropy density in hadron phase space as will be discussed in subsection 2.3. 
Not to be forgotten is the original observation about experimental data that has 
led to the data interpretation in terms of the sudden hadronization modet3: the 
strange quark fugacity as measured by emitted strange hadrons implies a source 
with freely moving strange s and antistrange s quarks such that (s — s) — 0, a point 
we will discuss in subsection 2.2. 

In subsection 3.3, we study in depth the phenomenon of strangeness enhance- 
ment and show that the rather precise analysis results we obtain are in excellent 
agreement with the theoretically computed strangeness yields, assuming formation 
of the QGP phase. We also show that, in the case of Pb-Pb, the explosive disintegra- 
tion of the dense QGP fireball leads to an overpopulation of the strangeness phase 
space abundance, and show that theoretical results are again in good agreement 
with the results of the data analysis. The initial temperatures for this agreement 
to occur are in the range 260 < T c h < 320, which values we obtain from models of 
collision dynamics. 

In section 4, we develop the theoretical method which leads to the finding of 
overpopulated strangeness phase space discussed in subsection 3.3. This, at a first 
sight surprising result, occurs due to early freeze-out of strangeness abundance in 
a rapid explosive evolution of the QGP fireball. We find a similar non-equilibrium 
result, in section 5, for RHIC condition in presence of transverse expansion which 
increases the speed at which QGP dilutes. We explore the dynamics of the phase 
space occupancy, rather than particle density, which allows us to eliminate much 
of the dependence on the dynamical flow effects by incorporating in the dynamics 
considered the hypothesis of entropy conserving matter flow and evolution. 

Exploiting the experience with SPS data analysis, we are able to consider, in 
section 5, strange particle production at RHIC. We obtain an unexpected parti- 
cle abundance pattern: during the hadronization of the baryon-poor RHIC-QGP 
phase there is considerable advantage for strangeness flavor to stick to baryons and 
antibaryons. This can be easily understood realizing that production of strange 
(anti)baryons is favored over production of kaons by the energy balance, i.e.: -B(A + 
7r) < £7(N+K). Moreover, at RHIC there are a high number of strange quarks per 
baryon available, and in this strangeness bath just a few (anti)baryons will manage 
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to emerge without strangeness content. We thus expect and find in detailed study in 
section 5 that hyperon production dominates baryon production, i.e., most baryons 
and antibaryons produced will be strange. We consider this result to be a unique 
consequence of the sudden QGP hadronization scenario observed at SPS and hope 
and expect that hyperon dominance should, when observed at RHIC, be generally 
accepted as proof of formation of the deconfined phase of nuclear matter. This 
phenomenon shows how much more pronounced will be the physics of strangeness 
in QGP at RHIC, as compared to the ten times lower SPS energy range. 

2. Contemporary Fermi Model of Hadron Production 
2.1. Phase space and parameters 

The relative number of final state hadronic particles freezing out from, e.g., a ther- 
mal quark-gluon source is obtained noting that the fugacity fi of the i-th emitted 
composite hadronic particle containing /c-components is derived from fugacities A& 
and phase space occupancies 7&: 

N^e- E ^h = e- E ^\{ lk X k . (1) 

We study chemical properties of light quarks u, d jointly, denoting these by a single 
index q and also consider chemical properties of strange quarks s . Thus as seen 
in Eq. (]]]), we study particle production in terms of five statistical parameters. In 
addition there is at least one matter flow velocity parameter. The six parameters 
which characterize the accessible phase space of hadronic particles made of light 
quarks 'q' and strange quarks 's' and their natural values, assuming a QGP source, 
are in turn: 

1) X s : The value of strange quark fugacity A s cab be obtained from the require- 
ment that strangeness balances, (N s — N s ) = , which for a source in which all s, s 
quarks are unbound and thus have symmetric phase space, implies A s = 1 . How- 
ever, the Coulomb distortion of the strange quark phase space plaxs an important 
role in the understanding of this constraint for Pb-Pb collisions ,E-a leading to the 
Coulomb-deformed value A s ~ 1.1 , as discussed in next subsection. 

2) 7 S : The strange quark phase space occupancy 7 S can be computed, Mild, will 
be studied in this review in detail within the framework of kinetic theory.li3lij For 
a rapidly expanding system the production processes will lead to an over-saturated 
phase space with 7 S > 1 . The difference between the two different types of chemi- 
cal parameters Xi and 7, is that the phase space occupancy factor 7^ regulates the 
number of pairs of flavor H\ and hence applies in the same manner to particles and 
antiparticles, while fugacity Xi applies only to particles, while A" 1 is the antiparticle 
fugacity. 

3) X q : The light quark fugacity X q , or equivalently, the baryochemical potential: 



/i B = 3T/lnA g , 



(2) 
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regulate the baryon density of the fireball and hadron freeze out. This density can 
vary dependent on the energy and size of colliding nuclei, and thus the value of X q is 
not easily predicted. However, since we know the energy per baryon content in the 
incoming nuclei, if we assume that the deposition of the baryon number and energy 
('stopping') is similar, we know the energy per baryon content in the fireballES This 
qualitative knowledge can be used in a study of equations of state applicable to the 
dense fireball to establish a constraint. 

4) 7 g : The equilibrium phase space occupancy of light quarks 7 g is expected to 
significantly exceed unity to accommodate the excess entropy content in the plasma 
phase. cjO There is an upper limit: 

lq < 7 * ee e-/ 2T , (3) 

which arises if all pions produced are simultaneously present, forming a Bose gas. 
We will address this effect in subsection 2.3. 

5) Tf. The freeze-out temperature Tf is expected to be not much different 
from the Hagedorn temperature Tji ~ 160MeVj23 which characterized particle 
production in proton-proton reactions. 

6) v c : The collective expansion velocity v c is expected to remain below the 
relativistic sound velocitylj: 

v c < 1/V3. (4) 

When the source emitting the free streaming particles is undergoing local collective 
flow motion, spectra of particles emitted are described by replacing the Boltzmann 
factor in Eq. (Q) by: 

7c = -T^i (5) 

a result which can be intuitively obtained by a Lorentz transformation between an 
observer on the surface of the fireball, and one at rest in the laboratory frame. 
Formal derivation of this and more elaborated results requires a considerably more 
precise frameworkEl 

The resulting yields of final state hadronic particles are most conveniently char- 
acterized taking the Laplace transform of the accessible phase space. This approach 
generates a function which in its mathematical properties is identical to the par- 
tition function. For example for the open strangeness sector we find, for the case 
v c = 0: 

VT 3 

\nZ s = — {{\ s \- x + X^X^^Fk + (X s X 2 q + A; 1 A- 2 ) 7s79 2 Fy 

+ (A 2 A 9 + a; 2 a- 1 ) 7s 2 7^h + (x 3 s + x; 3 ) 7 3 s F n } . (6) 

The integrated momentum phase space factors Fi for kaons i — K, single strange 
hyperons i — Y , doubly strange cascades i = 3 and triply strange omegas i — 
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are: 



Fi = J2 9 lj W(m lj /T) , W(x) = x 2 K 2 (x) . (7) 



gi j is the statistical degeneracy of each contributing hadron resonance 'j' of the 
kind V in the which comprise all known strange hadron resonances. K% is the 
modified Bessel function which arises from the relativistic phase space integral of 
the thermal particle distribution f(p) oc e -\/ m2 +P 2 /' r . It is important to keep in 
mind that: 

a) Eq. (^) docs not require formation of a phase comprising a gas of hadrons, but is 
not inconsistent with such a step in evolution of the fireball; in that sense it is not 
a partition function, but just a look-alike object arising from the Laplace transform 
of the accessible phase space, and 

b) the final particle abundances measured in an experiment are obtained after all 
unstable hadronic resonances 'j' are allowed to disintegrate, contributing to the 
yields of stable hadrons; 

c) the unnormalized particle multiplicities arising are obtained differentiating Eq. (^) 
with respect to particle fugacity. The relative particle yields are simply given by 
ratios of corresponding chemical factors, weighted with the size of the momentum 
phase space accepted by the experiment. For particles showing the^aame spectral 
shape comparison of normalization of m± spectra suffices, e.g., Ref. 113: 



2 (dss) 



A(dds) 

gi are the spin statistical factors of the states considered. Similarly: 

g-idis^K 2 



3 (dss) 



A(dds) 



(8) 



(9) 



When acceptance is limited to central rapidity, and significant flow is present con- 
siderable effort must be made to introduce appropriate phase space weights, 
d) In some experimental data it is important to distinguish the two light quark 
flavors as is in fact the case in the two above examples. This can be incorporated 
considering how the average light quark fugacity varies between both light quark 
speciesjld and assuming that the phase space occupancies are equal. 

We consider, for SPS energy range, the radial flow model, which is without doubt 
the simplest of the reasonable and expected matter flow cases possible, in view of the 
behavior of global observables seen in these experiments. As the results below show, 
this suffices to assess the impact of collective flow on the data analysis originally 
developed to be as little as possible sensitive to collective matter flow, even when 
particle yields in highly restricted regions of m±, y are considered. The collective 
source flow can completely change the shape of momentum distribution of particle 
produced, though of course it leaves unchanged the total particle yield, which is 
the integral sum of particle multiplicity over the entire phase space of the flow 
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spectrum. However, particles of different mass experience differing flow effects when 
mx, y acceptance cuts are present. Moreover, particles can freeze-out at slightly 
different conditions. In order to limit the influence of the practically unknown 
collective flow structure on particle yields in limited domains of the accessible phase 
space, we study compatible particle ratios: these are yield ratios obtained in a 
restricted domain of m± 1 y, for particles of similar mass and believed to have a 
similar interaction strength with the matter background. 

We now will address in turn two special topics, which slightly contradict expec- 
tations, and thus require more attention. Firstly, we review the the properties of 
the strange quark fugacity A s , which is sensitive to the possible asymmetry between 
strange and antistrange quarks in the source. The importance of this parameter is 
that it potentially helps distinguish the confined from deconfined phase: while in 
the baryon-rich confined phase the requirement of strangeness conservation implies 
that A s > f , in the deconfined phase the symmetry between phase space of strange 
and antistrange quarks implies A s ~ 1 . Following this, we address in more detail 
case of pions,E£l which is exceptional since we will be considering rather large val- 
ues of 7q > 1.5. As we shall see, the pion gas emerging from the QGP phase is 
strongly influenced by Bose correlation effects; in fact it is close to satisfying the 
Bose condensation condition. 



2.2. Coulomb force 



It has been recognized for a long time that the Coulomb force can be of con- 
siderable importance in the study of relativistic heavy ion collisions, ft plays an 
important role in the HBT analysis of the structure of the particle source.EJ We 
show that the analysis of chemical properties at freeze-out is also subject to this 
perturbing force, and in consideration of the precision reached in the study of par- 
ticle ratios, one has to keep this effect in mind. 

We consider a Fermi gas of strange and antistrange quarks, allowing that the 
Coulomb potential V established by the excess charge of the colliding nuclei dis- 
torts significantly the phase space. Within a relativistic Thomas- Fermi phase space 
occupancy model,t£l and allowing for finite temperature in QGP we haveEj: 



d 3 rd 3 p 
(2tt) 3 



1 



1 



1 



■iV(r))/T ! + 7s -l Ase (£( P ) + iV(r))/T 



(N s -N s ) =j g 

Rt 

(10) 

which clearly cannot vanish for V ^ in the limit A s — ► 1. In Eq. ( |To| ) the subscript 
Rf on the spatial integral reminds us that only the classically allowed region within 
the fireball is covered in the integration over the level density; E = yj m 2 + p 2 , and 
for a uniform charge distribution within a radius Rf of charge Zf\ 



V = 



3Z f e 2 
Z { e 2 



for r < Rf ; 
for r > Rf . 



(11) 
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One obtains a rather precise result for the range of parameters of interest to us 
using the Boltzmann approximation: 

(N. - N s ) = 7s \ [ g s -^- e ~ E / T X [ d?r \\ s e& - \?e~M . (12) 



' (2tt) 3 



Rt 



The Boltzmann limit allows us also to verify and confirm the signs: the Coulomb 
potential is negative for the negatively charged s-quarks with the magnitude of the 
charge, 1/3, made explicit in the potential terms in all expressions above. We thus 
haveO 

j" 7*6 ~ 

A s = A S A Q = 1 , A Q = f ^ . (13) 

•J R{ 

Aq < 1 expresses the Coulomb deformation of strange quark phase space. Aq is not 
a fugacity that can be adjusted to satisfy a chemical condition, since consideration of 
Aj, i = u, d, s, exhausts all available chemical balance conditions for the abundances 
of hadronic particles. The subscript Rf, in Eq. (|l3|), reminds us that the classically 
allowed region within the dense matter fireball is included in the integration over 
the level density. Choosing R f = 8fm, T = 140 MeV, m s = 200 MeV, noting that 
the value of 7 S is practically irrelevant as this factor cancels in Boltzmann approx- 
imation, see Eq. (|l2|), we find for Z f = 150 that the value A s = 1.10 corresponds 
to Rf — 7.9 fm. The Coulomb effect is thus relevant in central Pb-Pb interactions, 
while for S-Au/W/Pb reactions, similar analysis leads to a value A s = 1.01, little 
different from the value A s = 1 expected in the absence of the Coulomb phase space 
deformation. Another way to understand the varying importance of the Coulomb 
effect is to note that while the Coulomb potential acquires in the Pb-Pb case a 
magnitude comparable to the quark chemical potential, it remains small on this 
scale for S-Au/W/Pb reactions. 

2.3. Super-dense pion gas and chemical non- equilibrium 

For pions composed of a light quark-antiquark pair, the chemical fugacity is 7^, see 
Eq. (^). Thus the pion momentum space distribution has the Bose shape: 



ME) = _ 2 - - . E v = y/ml+p*. (14) 



l q 2 e 



The range of values for j q is bounded from above by the Bose singularity. When 
7 9 — > 7g, Eq. (^), the lowest energy state (in the continuum limit with p — > 0) will 
acquire macroscopic occupation and a pion condensate is formed. Formation of such 
a condensate 'consumes' energy without consuming entropy of the primordial high 
entropy QGP phase. On the other hand, as we shall see presently, when "/ q — > 7^ 
the entropy content of the pion gas initially grows! Thus while the development, 
directly from the QGP phase, of a pion condensate is not likely, the sudden had- 
ronization of entropy rich QGP should lead to the limiting value j q — > 7^, in order 
to more efficiently connect the entropy rich deconfined and the confined phases. 
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An interesting feature of such a mechanism of phase transition is that the chemical 
non-equilibrium reduces and potentially eliminates any discontinuity in the phase 
transition, which thus, in the experiment, will appear more like a phase transforma- 
tion without critical fluctuations, even if theory implies a 1st order phase transition 
for statistical equilibrium system. 

In Fig.[l], we show the physical properties of a pion gas as function of j q for a gas 
temperature T = 142 MeVE3 We see (solid line) that a large range of entropy density 
can be accommodated by varying the parameter j q . It is important to remember 
that in the hadronization of a quark-gluon phase it is relatively easy to accommodate 
energy density, simply by producing a few heavy hadrons. However, such particles 
being in fact non-relativistic at the temperature considered, are not effective carriers 
of pressure and entropy. However, as we see now in Fig.[l], the super-dense pion 
gas is just the missing element to allow a rapid hadronization process, since the 
entropy density is nearly twice as high at ^ q ~ 7^ than at j q = 1. Without this 
phenomenon one has to introduce a mechanism that allows the parameter VT 3 to 
grow, thus expanding either the volume V due to formation of so called mixed phase 
or invoking rise of T in so called reheating. 

The specific properties of the super-dense pion gas are shown Fig. ^|. In Fig. ^a, 
we relate the properties to the chemical equilibrium value j q = 1 and we also show 
that the Boltzmann approximation is not qualitatively wrong, as long as "f q < "f q ■ 
In Fig.J^b, we see the relative change in energy per pion, (inverse of) entropy per 
pion, and energy per entropy. Interestingly, we note that the entropy per pion drops 
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' 1 1 1 1 1 1 1 1 1 ' ' 1 1 1 1 1 1 1 1 1 1 ' 1 

0.5 1 1.5 0.5 1 1.5 

Fig. 2. Dependence of pion gas properties (jV-particle, i?-energy and S-entropy) density as 
function of jq for T = 142 MeV. a) ratios relative to equilibrium value 7 9 = 1; b) relative ratios, 
thus E/N, N/S and E/S. 



as 7 9 increases, and at the condensation point "f q = 7^, we can add pions without 
increase in entropy. We further note that a hadronizing gas will consume, at higher 
7 9 , less energy per particle, and that the energy per entropy is nearly constant. 

It is important to remember that if the hadronization process were adiabatic, 
allowing a full equilibrium relaxation, naturally 7, — > 1 would arise: as is implicitly 
well known the value 7 j — > 1 maximizes the entropv for a particle gas at fixed total 
energy, corresponding to the chemical equilibrium.!^ This result is easily found con- 
sidering Boltzmann pion gas, and recalling in some detail the definition of entropy 
since the standard equilibrium expressions do not apply: 

Sb,f = J^0[±(l±f(x,p))ln(l±f(x,p))-f(x,p)\nf(x,p)}, (15) 

^ V J(^k fip)He/mh (16) 

where aside of 'B,F' (Bose, Fermi) also the Boltzmann limit for a homogeneous 
spatial distribution is shown explicitly. 

Evaluating in Boltzmann limit the particle number and energy, we find that the 
factor 7 2 becomes a normalization factor which describes the average occupancy of 
the phase space relative to the equilibrium value, and for entropy, we also find a 
logarithmic term: 

N = 7g 2 iV| cq -> aV^T 3 , (17) 
E = 7g 2 £| oq -> 3aV 1 2 q T i , (18) 
S = j 2 q S\ eq + In ( 7 " a ) 7 2 7V| oq -> iaV^T 3 + In ( 7 " a ) aV-fiT 3 . (19) 
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For massless pions, a = g/ir 2 , with pion degeneracy g — 3 . Setting E = Const., we 
eliminate T and find that the entropy as function of j q varies according to: 

%=Const. «7l(4-ln 79 2 ), (20) 

which has a very weak maximum at j q = 1; note that at 7 = 1.4, the entropy is at 
98.3% of the value at 7 = 1. At this point it is important to realize that the chemical 
equilibrium is much better defined for the more familiar case of a fixed temperature 
bath (and not an isolated fixed energy fireball discussed above): consider the free 
energy T of the (non-interacting) relativistic gas at fixed temperature T . Since 
T = E — TS, we combine Eqs. (|17|-|Ig|) and obtain in the Boltzmann limit, 

^ = - a yT 4 7 9 2 [l + ln( 7 - 2 )] , (21) 

which has a minimum for chemical equilibrium value j q = 1. However, one now 
finds that a change by factor 1.4 in j q) at fixed T leads to a change by 35% in the 
value of the free energy and even a greater change in entropy. 

We conclude that, in adiabatic condition, the fireball would evolve to the max- 
imum entropy equilibrium case 7; = 1,2 = q, s, the gain in entropy for an isolated 
system is in this limit very minimal. Thus for a system with rapidly evolving vol- 
ume, we will in general find more effective paths to increase entropy, other than 
the establishment of the absolute chemical equilibrium. Hence the values we report 
7, 7^ 1,2 = q, s, are consistent with the present day understanding of explosive 
evolution of the hadronic matter fireball. 

In a systematic study of the relevance of different physical parameters, the chem- 
ical non-equilibrium at hadron freeze-out has been shown to be a required ingredient 
in order to arrive at a precise interpretation of the experimental results on parti- 
cle ratios ffl obtained at CERN. This is best seen considering, the results for the 
statistical parameters obtained for the S-Au/W/Pb collisions,^ and the associated 
total statistical error, 

2 = J2j( R lh ~ -ftcxp) 2 , 22 < 

XT ~ (ARL p ) 2 ' 
which are presented in table [|. We clearly see the gain in physical significance that 
is accomplished as chemical non-equilibrium is allowed for by releasing the fixed 
value 7j = 1, first for strange and next, light quarks. We also observe that allowing 
for X s 1 does not lead to an improvement in statistical significance, since the data 
is compatible with this value expected for the deconfined QGP. Similar systematic 
study has also been completed for the Pb-Pb system,CJ reconfirming the need to 
use 7^ 7^ 1 in the data analysis. 

The errors in the results shown in table and in results that follow below, 
are one standard deviation errors arising from the propagation of the experimental 
measurement error. However, these errors are meaningful only when the theoretical 
model describes the data well, as is the case for last entry line in table [j] when we 
allow light quark chemical nonequilibrium, 7 9 7^ 1 . 
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Table 1. Statistical parameters obtained from fits of S-Au/W/Pb data without enforcing strangen- 
ess conservation, and not considering flow effect — only compatible particle ratios were considered. 
Asterisk (*) means a fixed (input) value. See text for an explanation of result errors. 



T f [MeV] 


\ 




Is 


7 9 


X^/dof 


145 ± 3 


1.52 ± 0.02 


1* 


1* 


1* 


17 


144 ± 2 


1.52 ± 0.02 


0.97 ± 0.02 


1* 


1* 


18 


147 ± 2 


1.48 ± 0.02 


1.01 ± 0.02 


0.62 ± 0.02 


1* 


2.4 


144 ± 3 


1.49 ± 0.02 


1.00 ± 0.02 


0.73 ± 0.02 


1.22 ± 0.06 


0.90 



3. Strange Hadron Data Analysis 
3.1. Particle yields 

The available compatible particle yield ratios (excluding presently Vl and fl, see 
subsection 3.4) are listed in table ||, top section from the experiment WA97, for 
p± > 0.7 GeV, within a narrow Ay = 0.5 central rapidity window. Further below 
are shown results from the large acceptance experiment NA49, extrapolated by the 
collaboration to full 4tt phase space coverage. We first fit 11 experimental results 
shown in table and then turn to include also the m^-slope into this consideration, 
and thus have 12 data points. The total error \\ f° r the four result columns is shown 
at the bottom of this table along with the number of data points W, parameters 'p' 
used and (algebraic) redundancies V connecting the experimental results. For r =/= 
it is more appropriate to quote the total Xti with a initial qualitative statistical 
relevance condition x^/(N — p) < 1. 

The first theoretical columns refer to results without collective velocity v c (sub- 
script 0) the three other were allowing for v c (subscript v c ). In column three, 
superscript 'sb' means that A s is fixed by strangeness balance and, in column four, 
superscript 'sc' means that 7 9 = 7^ = e m ^l 2T s , that is 7 g is fixed by its upper limit, 
the pion condensation point. AlLresults shown account for slightly higher value of 
the ratio h~ /B recently reportedES; B is here the number of baryon participants and 
h~ = tt~ + K~ + p is the yield of stable negative hadrons comprising as indicated 
pions, kaons and antiprotons. 

First we note that all columns in table |^ represent physically acceptable result 
for the Pb-Pb collision system: 

a) presence of collective flow (three last columns) leads to very similar compatible 
particle ratios, even though improvement of xt occurs when v c 7^ is allowed for; 

b) the highest confidence result is found just when the light quark phase space 
occupancy assumes a value at below the pion condensation point; 

c) strangeness conservation (enforced in second last column) is naturally present, 
enforcing it does not change in any way the results for particle multiplicities. 

Allowing radial flow not only improves the capability to describe the data, but 
it allows us to study m± particle spectra, which offer another independent measure 
of flow, and confirm the value of v c — when considering v c along with Tj_, the 
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Table 2. WA97 (top) and NA49 (bottom) Pb-Pb 158A GeV particle ratios compared with theo- 
retical results, see text for explanation. 



Ratios 


I 


Le 


f. 


Exp. Data 


Pb| 


Pb|„ 


Pb|f 


Pb|* 


S/A 




34 




0.099 ± 0.008 


0.104 


0.103 


0.105 


0.103 


s/A 




34 




0.203 ± 0.024 


0.214 


0.208 


0.209 


0.206 


A/ A 




34 




0.124 ± 0.013 


0.124 


0.125 


0.124 


0.125 






34 




0.255 ± 0.025 


0.256 


0.252 


0.248 


0.251 


(S+3) 
(A+A) 




36 




0.13 ± 0.03 


0.126 


0.122 


0.124 


0.122 






35 




11.9 ± 1.5 


14.2 


13.3 


13.0 


13.4 


K+/K- 




33 




1.80± 0.10 


1.80 


1.82 


1.78 


1.83 


p/P 




BI 




18.1 ±4. 


17.3 


16.7 


16.6 


16.6 


A/p 




37 




3. ± 1. 


2.68 


2.11 


2.11 


2.11 


K° s /B 




SS 




0.183 ± 0.027 


0.181 


0.181 


0.163 


0.188 


h~/B 




36 




1.97 ± 0.1 


1.96 


1.97 


1.97 


1.96 




At 




3.6 


2.5 


3.2 


2.6 




N;p;r 




11;5;2 


12;6;2 


12;5;2 


12;5;2 



Table 3. Expcrime«*»J and theoretical spectra inswse slopes T t ^. Left Pb-Pb results from 
experiment WA97E3H; right S-W results from WA85.E3 





T| b [McV] 


T&> [McV] 


[MeV] 


[MeV] 


rpK U 


230 ± 2 


241 


219 ± 5 


215 


rpA 


289 ± 3 


280 


233 ± 3 


236 


rpA 


287 ± 4 


280 


232 ± 7 


236 


T 3 


286 ± 9 


298 


244 ± 12 


246 




284 ± 17 


298 


238 ± 16 


246 



inverse slope of the m± spectra, we have one parameter and several spectral inverse 
slopes of particles considered. However, we will in the first instance assume that we 
have just one additional data point and we proceeded as follows: for a given pair 
of values Tf and v c we evaluate the resulting m± particle spectrum and analyze it 
using the spectral shape and kinematic cuts employed by the experimental groups. 
Once we find values of Tf and v c , we study again the inverse slopes of individual 
particle spectra and obtain an acceptable agreement with the experimental T| as 
shown in left section of table ||. We have considered in the same framework the 
S-induced reactions, and the right^ection of table || shows also a good agreement 
with the WA85 experimental data.E3 

We have updated the experimental Pb-Pb results shown in table || with the cur- 
rent high precision results.Ell However, the theoretical results shown were obtained 
earlier for slightly different results with larger error bars, and we hope to reevaluate 
these results in the near future. To model these slopes theoretically, one needs to 
remember that the vast majority of A and A is a decay product of S° and S°, A* 
and A* and S and S. Consequently given the precision of the (inverse) slopes pre- 
sented, in order to model the A and A spectra one will need to consider the effect 
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of hadron cascading, which introduces uncertainty arising from a dependence on 
unmeasured yields. However, given the current availability of quite precise S and S 
slopes, and the fact that these particles are rarely decay products of other hadronic 
resonances, we will in future use these slopes as the spectral data point input in the 
data analysis studies. As result, we anticipate a slight reduction in the collective 
velocity within the errors shown below. 

3.2. Chemical freeze-out properties 

The six parameters (T/, v c , X q , X s , j q , j s ) describing the particle abundances are 
shown in the top section of table || We also show in the last column the best result 
for S-induced reactions, where the target has been W/Au/Pb.o All results shown 
in table [| have convincing statistical confidence level. For the S-induced reactions 
the number of redundancies r shown in the heading of the table ^|is large, since the 
same data comprising different kinematic cuts has been included in the analysis. 

Within error, the freeze-out temperature Tf ~ 143 ± 3 MeV, seen in table f|, 
is the same for both the S- and Pb- induced reactions, even though the chemical 
phase space occupancies differ greatly. Such a behavior is expected in view of 
the similarity of the energy content in the collision in both reaction systems, but 
greatly differing collision geometry. We find that the variation in the shape of 
particle mj_-spectra is fully explained by a change in the collective velocity, which 
rises from = 0.49 ± 0.02 to v^ h = 0.54 ± 0.04 ~ 1 / V3 = 0.577. The value of light 
quark fugacity X q implies that baryochemical potential is /i^ b = 203 ± 5 > n% = 
178 ±5 MeV. As in S-induced reactions where A s = 1, now in Pb-induced reactions, 
a value A^ b ~ 1.1 characteristic for a source of freely movable strange quarks with 
balancing strangeness, i.e., A s = 1, is obtained, see Eq. ([l3|). 

Further evidence for low chemical freeze-out temperature is contained in the 
m^-particle spectra considered in subsection 3.1. Our approach offers a natural 
understanding of the equality of the m^-slopes of the strange baryons and an- 
tibaryons considered which arises because within the sudden hadronization model 
both these particles emerge free-streaming from QGP. In the hadron based micro- 
scopic simulations this behavior of m^-slopes of baryons and antibaryons arises 
from fine-tuning of the particle-dependent freeze-out times.fEl On the other hand, in 
such a microscopic study one finds in view of the small reaction cross sections that 
f2 and Q could freeze out somewhat sooner than the other hadrons, and thus ps/puld 
have a softer spectrum as also confirmed in direct hadronization simulationsJla We 
will return to this point just below. The reader should keep in mind that since we 
find a rather low chemical freeze-out temperature, and can explain the m± spectra 
well based on this value, the implied kinetic (collision) freeze-out temperature must 
be rather similar to the chemical freeze-out. 

In the bottom section of table |4|, we also see the energy and entropy content 
per baryon. The energy per baryon seen in the emitted hadrons is nearly equal to 
the available specific energy of the collision (8.6 GeV for Pb-Pb, 8.8-9 GeV for S- 
Au/W/Pb). This implies that the fraction of energy deposited in the central fireball 



Observing Quark- Gluon Plasma with Strange Hadrons. 19 



Table 4. In heading, we present the total quadratic relative error x^, number of data points 
N, parameters p and redundancies r; in the upper section: statistical model parameters which 
best describe the .experimental results for Pb— Pb data, and in last column for S-Au/W/Pb data 
presented in Ref.|fj. Bottom section: specific energy, entropy, anti-strangeness, net strangeness 
of the full hadron phase space characterized by these statistical parameters. In column two, we 
fix A s by requirement of strangeness conservation, and in column three, we fix 7 9 at the pion 
condensation point j q = jS. 





Pb\ v 


Pb|f 


Pb|~ 


su 


Xt; N;p;r 


2.5; 12; 6; 2 


3.2; 12; 5; 2 


2.6; 12; 5; 2 


6.2; 16;6;6 


T f [MeV] 


142 ± 3 


144 ± 2 


142 ± 2 


144 ± 2 


v c 


0.54 ± 0.04 


0.54 ± 0.025 


0.54 ± 0.025 


0.49 ± 0.02 


\ 


1.61 ± 0.02 


1.605 ± 0.025 


1.615 ± 0.025 


1.51 ± 0.02 




1.09 ± 0.02 


1.10* 


1.09 ± 0.02 


1.00 ± 0.02 


lq 


1.7 ± 0.5 


1.8± 0.2 


7 c* = e m,/2T, 


1.41 ± 0.08 




0.79 ± 0.05 


0.80 ± 0.05 


0.79 ± 0.05 


0.69 ± 0.03 


E f /B 


7.8 ± 0.5 


7.7 ± 0.5 


7.8 ± 0.5 


8.2 ± 0.5 


S f /B 


42 ± 3 


41 ± 3 


43 ± 3 


44 ± 3 


s f /B 


0.69 ± 0.04 


0.67 ± 0.05 


0.70 ± 0.05 


0.73 ± 0.05 


(*f - s f )/B 


0.03 ± 0.04 


0* 


0.04 ± 0.05 


0.17 ± 0.05 



must be nearly (within 10%) the same as the fraction of baryon number. The small 
reduction of the specific entropy in Pb-Pb compared to the lighter S-Au/W/Pb 
system maybe driven by the greater baryon stopping in the larger system, also seen 
in the smaller energy per baryon content. Both collision systems freeze out at the 
same energy per unit of entropy, 

E/S = 0.185 GeV. 

There is a loose relation of this universality in the chemical freeze-out condition 
with the suggestion made recently That particle freeze-out occurs at a fixed energy 
per baryon for all physical systems,0 considering that the entropy content is related 
to particle multiplicity. The overall high specific entropy content we find_agrees well 
with the entropy content evaluation we made earlier for the S-Pb case.c3 The high 
entropy content is observed in the final hadron state in terms of enhanced pion yield. 
Thus the ratio of K + /ir + is combines these two effects and is not a good indicator 
of new physics, even though this relatively simple observable continues to attract 
attention.113 It would have been more useful if systematic studies of strangeness 
production and enhancement were to offer as result of their analysis the strangeness 
yields per participating baryon number. 

The large values of "f q > 1, seen in table ||, imply as discussed earlier that there is 
phase space over-abundance of light quarks, which receives contribution from, e.g., 
gluon fragmentation at QGP breakup. j q assumes in the data analysis a value near 
to where pions could begin to condense£3 Eq. (||). This result is consistent with 
the expectations for hadronization of an entropy rich quark gluon plasma, as we 
discussed above in subsection 2.3. We found by studying the ratio h~ / B separately 
from other experimental results that the value of 7 9 ~ 7^ is fixed consistently and 
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independently both, by the negative hadron (h~), and the strange hadron yields. 
The unphysical range j g > ~ 1.63 can arise (see column Pb|i)f) since, up to this 
point, we had used only a first quantum (Bose/Fermi) correction. However, when 
Bose distribution for pions is implemented, which requires the constraint "f q < 7^, 
we obtain practically the same results, as shown in the third column of table |l|. 

3.3. Strangeness enhancement 

We show, in the bottom section of table |], the specific strangeness content, Sf/B 
along with specific strangeness asymmetry (s / — s /) / B seen in the hadronic particles 
emitted. In the data analysis the requirement that the number of s and s quarks 
in hadrons is equal is in general not enforced. We see that in lower portion of 
table |^ that this result is found automatically for the symmetric Pb-Pb collision 
system. However, a 3.5 s.d. effect is seen in the asymmetrical S-Au/W/Pb system 
Though the errors which we derive from the experimental data are small, there could 
be in this asymmetric system a considerable systematic experimental error due to 
data extrapolations made in presence of a significant spectator matter component, 
coupled with theoretical error from the varying CM-rapidity. On the other hand, 
the consistency of the Pb-Pb and S-Au/W/Pb results suggest that this asymmetry 
is possibly a real effect, thus the unseen balance of strangeness could be hidden 
in a residual (strange) quark matter nugget, which is escaping detection. Such 
strangeletts could in principle form, since in the hadronization of the S-Au/W/Pb 
deconfined sy s t em th e hadron phase space is asymmetric, which leads to strangeness 
distillation.O0@E 

One of the important quantitative results of this analysis is shown in the bottom 
section of table ||: the high yield of strangeness per baryon, Sf/B ~ 0.7. We now 
proceed to verify if this yield is in agreement with the pedictions made over the 
past 20 years. Perhaps more by chance than design, this analysis result is in agree- 
ment .Bath the first calculations of strangeness production employing perturbative 
QCD,li3 where the value N s /B = n s jv — 0.7 is reached for the plasma temperature 
of 300 MeV as shewn there in Fig. 3. Since, considerably more refined methods have 
been developedJiS and these are in excellent agreement with results of the analysis of 
experimental results. In view of the high precision reached in this data analysis, we 
have recomputed the theoretical yield taking for the QCD parameters values gen- 
erally accepted today: a s (Mz) = 0.118 and m s (lGeV) = 200 MeV, correpsonding 
to m s (M z ) = 90 MeV. 

In table g, we summarize for three collision systems we consider S-Au/W/Pb, 
Ag-Ag, Pb-Pb the key input parameters used in computing the result for N s /B 
shown below in Fig.^. The first entry line gives the central collision particpant 
numbers for the three systems considered. Next, in table ^, we see the initial tem- 
perature T c h which the evaluation of strangeness production requires as input. T c ^ 
is the temperature at the time when light quarks and gluons reach equilibrium. To 
obtain this value, we compute the collisional pressure and set it equal to thermal 
pressure at the time the fireball begins to expand.EffiJ To do this we need the (mo- 
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Table 5: Input to the strangeness production computation in QGP. 





S-Au/W/Pb 


Ag-Ag 


Pb-Pb 


participants 


B 


90 


180 


360 


q, G equilibration 


T ch [MeV] 


260 


280 


320 


stopping 


V 


52% 


54% 


57% 


scale expansion velocity 


V c 


0.49c 


0.52c 


0.54c 



mentum. energy) stopping fractions r\ here taken from NA35/NA49 experimental 
results,E3 (except for interpolation for Ag-Ag, the dotted line in Fig.||). The last 
line in table |^ addresses the expansion dynamics we use: we employ the observed 
freeze-out expansion velocity v c as given in the top section of table ^. We assume 
that each local volume expands its size scale R at this local velocity, and we consider 
the process to be entropy conserving, hence we use 

R 3 T 3 = Const. to obtain the 

time dependence of local fireball temperature. 

We obbtain the result for N s /B shown in Fig. [| as function of the specific energy 
available in the fireball E/B, for the three collision systems S-Au/W/Pb (short- 
dashed line), Ag-Ag (long dashed) and Pb-Pb (solid line). Since we compute the 
intial temperature from the collision energy our approach allows us to extrapolate 
as function of E/B, assuming that the stopping fraction for the collisional pressure 
is known. When we keep the stopping fraction constant and as given at the 160- 
200A GeV collision energy, we find the results shown in Fig.|3[ However, a constant 
stopping underestimates the intial temperature at lower collision energy, where we 
would expect higher stopping, and it overestimates the initial temperature at higher 
collision energies, where we would expect smaller stopping, thus we believe that the 
slope of the result we present in Fig.|| is too steep. We will be able to improve on 
this result after the behavior of stopping as function of collision energy has been 
understood. 

In Fig.|^, the solid square is the result of the analysis for S-Au/W/Pb system, 
and open square for Pb-Pb as shown in the lower section of table |]. We note that 
the reason that the available energy E/B in the fireball is the dominant parameter 
controlling strangeness yield is the cancellation of effect of higher initial temperature 
in the larger, more stopping systems, by the faster explosion of such a system, which 
leaves less time for strangeness production. We note that even though we did not 
analyze here the S-S system, for which case we would need to adapt the method 
to allow significant longitudonal flow, it is understood that the available fireball 
energy and strangeness content per baryon is higher in S-S 200 A GeV interactions,^ 
consistent with the results shown in Fig.||. 

This high strangeness yield corresponds to (above) equilibrium abundance phase 
space occupancy in hadronization. In the top section of table Q the ratio Js/jq — 
0.8, which corresponds (approximately) to the parameter j s when j q = 1 has been 
assumed. We observe that 7^ b > 1. This strangeness over-saturation effect could 
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E/B [GeV] 

Fig. 3. QGP fireball specific per baryon strangeness abundance as function of E/B energy per 
baryon in the fireball, for running a s (Mz) = 0.118 and m B (lGeV) = 200 MeV, correpsonding to 
m s {Mz) = 90 MeV. Solid lines: for Pb— Pb, stopping 57% and 360 particpants, long-dashed lines: 
for Ag— Ag, stopping 54% and 180 participants, short-dashed lines: for S— Au/W/Pb stopping 52% 
and 90 participants in QGP fireball. The solid square is the result of an analysis for S-Au/W/Pb 
200A GeV reaction system, and open square for Pb— Pb at 158A GeV, as shown in the lower section 
of table ^ 



arise from the effect of gluon fragmentation combined with early chemical equili- 
bration in the QGP, 7 s (t < tj) ~ 1. The ensuing rapid expansion preserves this 
high strangeness yield, and thus we find the result 7 S > 1 , as we reported in Ref. 
This high phase space occupancy is one of the requirements for the enhancement 
of multi-strange (aiiti)baryon production, which is an important hadronic signal of 
QGP phenomena. til 

We compare this result of data analysis, in quantitative manner, with the the- 
oretical computation of 7 S which is easily obtained from the above study of total 
strangeness production, as we only need to divide the total momentary strangen- 
ess yield by the expected equilibrium abundance, for which we choose to consider 
ideal gas of strange quarks with QCD running mass m s (fi = 5.5T) . The factor 5.5 
converts the value of T into the appropriate scale /j, of energy at which the kinetic 
equilibrium distribution is formed, and we note that m s (T — 182 MeV) = 200 MeV. 
The effect of QCD running influences the agreement between theory and experiment 
at the level of 10-15%. The result is shown in figure ||, right as function of time 
t for the 160-200A GeV collision systems and left as function of temperature T. 
Horizontal dotted line refers to equilibrium phase space occupancy, and the vertical 
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Fig. 4. Evolution of strangeness phase space occupancy f a , a) as function of time t (on left) and 
b) as function of temperature T (on right) for the 160— 200A GeV collision systems. For further 
explanation see legend to figure H and text. 



line indicates expected freeze-out condition at Tf — 143 MeV. 

Solid dots in figure ||a) show where this freeze-out temperature occurs as function 
of time t. The analysis point uncertainty in freeze-out time is obtained assuming 
that, in isentropic evolution, size scale R and temperature satisfy RT — Const., 
and thus with v c = dR/dt\f we find: 

M = -^=AT. (23) 

The chemical freeze-out occurs at about lOfm/c and 13fm/c, after onset of the 
collision, allowing for about 1 fm/c initial time r c h for the two systems S-Au/W/Pb 
and Pb-Pb respectively. Considering that the expansion velocity has been 0.5c and 
0.64c respectively we obtain an estimate of the chemical freeze-out radius R^ ~ 5fm 
and R^ h ~ 8fm, the latter value is in excellent agreement with the discussion of the 
Coulomb effect presented in section 2.2. 



3.4. The Omega riddle 

The QGP formation and sudden disintegration model we have described above 
has natural limitations. When we attempt to describe within this approach most 
rarely produced particles, there is the potential for under-prediction of experimen- 
tal results, which could receive contributions from other more effective production 
mechanisms. In this context, the most rarely produced hadron is the triply strange 
fi(sss) and fl(sss) which are the heaviest stable hadrons, Mq = 1672 MeV. The 
phase space for f2 is more than 10 times smaller than that for S at the conditions 
of chemical freeze-out we have obtained, and f2 production pattern can thus be 
altered by processes not implemented in the one stage fireball model used to analyze 
the data. 
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Table 6: Pb-Pb 158A GeV particle ratios involving f2 and f2, compared to theo- 
retical expectations for Tf = 143 MeV. * reminds us that statistical parameters are 
fixed by other particle yields. Last column presents results allowing a 11% H and 
E shadow, see text for details. 



Ratios Re 


\ Exp. Data 


v = 


V c 


11 % 


n/E 


m 


0.192 ± 0.024 


0.078* 


0.077* 


0.186 


Q/E 




0.27 ± 0.06 


0.17* 


0.18* 


0.28 


Q/Q 




0.38 ± 0.10 


0.57* 


0.60* 


0.38 


(n+Ti) 

(H+E) 


E 


0.20 ± 0.03 


0.10* 


0.10* 


0.206 



When we attempted to describe along with the other hadrons the yields of £1 
and fl within the single stage freeztput model, we indeed have discovered consid- 
erable loss of physical significance. CJ Already for the S-Au/W/Pb case, we have 
found that a more reliable description of the data arises if we did not consider the 
qualitative Omega yields available.E-3 For the parameters as reported, we find in 
the Pb-Pb reactions that we under-predict the f2 and f2 yields by about factor 2. 
The experimental results are shown in the first three columns of table ^ and the 
theoretical yield computed using the Fermi-2000 model with parameters fixed by 
other particle abundances are shown in columns 4 and 5: we see that the presence 
of radial flow (v = v c ) has a minimal impact on the relative yields, compared to 
the case without radial flow (v = 0). To put this result into proper perspective, 
consider that we find within the sudden hadronization of QGP with uncorrelated 
strange quarks in the deconfincd phase an enhancement of Q and yields 'only' 
by factor 10 as compared to what is expected from extrapolation of p-A reactions. 
However, the experiment reports an enhancement by factor 15-20. Such a 'failure' 
is in fact confirming the early expectations that f2 and f2 yields are the best si 
nature of deconfinement, considering the possibility of strange quark clustering Jl 
In fact it is a bit surprising how well this early prediction works, and this requires 
further study to understand more precisely what exactly this means. 

Several groups have noted, studying the microscopic evolution of and f2, that 
due to low reaction cross section th£v_decouple from hadron background somewhat 
sooner than all the other hadrons .EallJ An early chemical freeze-out would impact 
statistical yields of and fl greatly. To augment the f2 and il yields by factor k, it is 
sufficient to take an incrementally ST higher freeze-out temperature, as determined 
from study of the phase space: 

In k , 

ST ~ T — . (24) 

M n /T 1 ; 

Thus in order to increase the yields by a factor 2 the Q and fl freeze-out would need 
to occur at Tq = 150 MeV rather than at Tf — 143 MeV. Since the temperature 
drops as the explosion of the fireball develops, this higher freeze-out temperature 
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means an earlier in time freeze-out. 

Even if the required staging in time of hadron production is apparently small, 
a consistent picture requires fine-tuning and it seems unnatural, considering that 
all the other particles are perfectly consistent with just one sudden freeze-out con- 
dition. Pursuing other alternatives, we note that fl and enhancement is caused 
by strangeness pre-clustering in the deconfined phase which would enhance mul- 
tistrange hadrons, but most prominently and noticeable enhance the phase space 
suppressed SI and SI. In this context, it is interesting to note that the missing yield 
is not symmetric: as seen in table [| we miss in relative terms more SI than SI . In- 
terestingly, the missing yield is exactly proportional to the yield of 5 and H and the 
best description of all particle yields, including all SI and SI is arrived at describing 
what is missing as proportional (11 %) to the 2 and 2 yield, this is shown in the 
last column '11 %' of table ^. It is now easy to propose a model that would lead 
just to this result: there are colored di-strange quarks clusters at hadronization and 
when their color strings break 5 and SI are produced. This imprints a 'shadow' of 2 
and 2 in the SI and Sl-abundance. While this works for f2 and SI, we find that this 
mechanism is not compatible with the other particle abundances, in other words a 
similar 'shadow' of A and A in the 2 and 2 channel seems unacceptable. Thus this 
mechanism would work only if pairing of strange quarks would be significant near 
to phase transition. Current models of 'color super conductivity' support such a 
clustering mechanism for additional SI and SI enhancement, though detailed studies 
are still in progress. 113 

We have also explored the possibility that unknown SI* and SI* resonances con- 
tribute to the SI and fl yield, but we were not able to find a good set of pa- 
rameters for these hypothetical resonances. Moreover, this hypothesis implies a 
baryon-antibaryon symmetric contribution in the sense that both SI and SI yields 
are multiplied by the same factor. However, the missing yield is clearly also baryon- 
antibaryon asymmetric — thus despite several ad-hoc parameters the model descrip- 
tion remains poor. 

We note that earlier statistical descriptions of SI and SI yields have not been 
sensitive to the problems we described.oO In fact as long as the parameter 7 g is 
not considered, it is not possible to describe the experimental data at the level of 
precision that would allow recognition of the S7 and SI yield as a problem for the 
statistical Fermi phase space model. 

4. Kinetic Strangeness Production 

In so me compu tational details, the methods to describe strangeness production 
differ .EjOSH This leads to different expectations regarding chemical equilibra- 
tion of quark flavor at RHIC energies, with some authors finding marginal at best 
chemical equilibration. We therefore develop in more detail the computational ap- 
proach which is consistent with the SPS-energy scale results discussed in previous 
sections. One important difference to the earlier work is that the two loop level 
running of QCD parameters for both coupling strength a s and strange quark mass 
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m s is used. olm z = 0.118 is assumed as determined at the p. = Mz° energy scale. 
Another improvement is that an entropy conserving explosive flow of matter is in- 
corporated directly into the dynamical equations describing the evolution of stran- 
geness phase space occupancy. This approach is entailing significant cancellations 
in the dynamical equations and the only model dependence on matter flow which 
remains is the relationship between the local temperature and local proper time. In 
consequence, a relatively simple and physically transparent model for the evolution 
of the phase space occupancy 7 S of strange quarks in the expanding QGP can be 
studied. 

We use two assumptions of relevance for the results we obtain: 

• the kinetic (momentum distribution) equilibrium is reached faster than the chem- 
ical (abundance) equilibriurrO'E3; 

• gluons equilibrate chemically significantly faster than strangeness. lB 

The first assumption allows us to study only the chemical abundances, rather than 
the full momentum distribution, which simplifies greatly the structure of the master 
equations; the second assumption allows us to consider the evolution of the stran- 
geness population only after an initial time r c h period has passed: r c h is the time 
required for the development to near chemical equilibrium of the gluon population, 
and the corresponding temperature T c h is the initial condition we need to compute 
the evolution of strangeness. Aside of T c h, the strange quark mass m s introduces 
the greatest uncertainty that enters strangeness yield calculations based on resumed 
perturbative QCD ratesO The overpopulation of the strangeness phase space, seen 
before in section 3 in SPS data, arises in particular for T c h > 250 MeV and values 
of strange quark mass m s (lGeV)~ 200 ± 20 MeV. 

In view of these assumptions the phase space distribution f s can be characterized 
by a local temperature T(x, t) of a (Boltzmann) equilibrium distribution , with 
normalization set by a phase space occupancy factor: 

fM^t))^is(T)fr(p;T). (25) 

Eq. ( p5| ) invokes in the momentum independence of 7 S the first assumption. More 
generally, the factor 7^, i — g,q,s,c, allows the local density of gluons, light quarks, 
strange quarks and charmed quarks, respectively to evolve independently of the local 
momentum shape. With variables (t, x) referring to an observer in the laboratory 
frame, the chemical evolution can be described by the strange quark current non- 
conservation arising from strange quark pair production described by a Boltzmann 
collision term: 



ll= d£a dVp. _ 1 2fl , )gg - 

p q (t) Pq (t)(av)f- sS - Ps (t) Ps (t) (av) s T ^ a3 >™. (26) 



The factor 1/2 avoids double counting of gluon pairs. The implicit sums over 
spin, color and any other discreet quantum numbers are combined in the particle 
density p = ^ s J d 3 p /, and we have also introduced the momentum averaged 
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production/annihilation thermal reactivities (also called 'rate coefficients'): 
, . I d 3 Pi J d 3 p 2 a 12 v 12 f(pi,T)f(p 2 , T) 

{aVrcl)T = /^/«(?i,T)/fer) ■ (27) 

f(j>i,T) are the relativistic Boltzmann/Jiittner distributions of two colliding parti- 
cles of momentum pi, i = 1,2. 

The current conservation used above in the laboratory 'Eulerian' formulation can 
also be written with reference to the individual particle dynamics in the so called 
'Lagrangian' description: consider p s as the inverse of the small volume available 
to each particle. Such a volume is defined in the local frame of reference for which 
the local flow vector vanishes v(x, t)|i OC al = 0. The considered volume SVi being 
occupied by small number of particles SN (e.g., SN = 1), we have: 

SN S = Ps 5Vi ■ (28) 

The left hand side (LHS) of Eq. (E6I) can be now written as: 



dp s dvps _ 1 dSN s _ dp s 1 dSVi 

dt dx ~ SVi dt ~ dt ps 5Vi dt ' ( ' 

Since SN and SVidt are L(orentz)-invariant, the actual choice of the frame of refer- 
ence in which the right hand side (RHS) of Eq. ( |29| ) is studied is irrelevant and we 
drop henceforth the subscript I. 

We can further adapt Eq. ( p9| ) to the dynamics we pursue: we introduce p£°(T) 
as the (local) chemical equilibrium abundance of strange quarks, thus p = 7 s /c^°. We 
evaluate the equilibrium abundance SN^° = 5Vpf(T) integrating the Boltzmann 
distribution: 

SN™ = [SVT 3 ]^z 2 K 2 (z) 1 z=^. (30) 

We will below use: d\z v ' K v (z)]/dz — —z v K v ^\, where K v is the modified Bessel 
function of order v. The first factor on the RHS in Eq. ( |30| ) is a constant in time 
should the evolution of matter after the initial pre-thermal time period To be entropy 
conserving,!! and thus 5VT 3 = SVqTq = Const. . We now substitute in Eq. ( |29| ) 
and obtain 

d Ps , dv Ps _ fpT (dj 1 + , (31) 



dt dx \dT T K 2 {z) j 

where T = dT/dt. Note that, in Eq. (|3l|), only a part of the usual flow-dilution 
term is left, since we implemented the adiabatic volume expansion, and study the 
evolution of the phase space occupancy in lieu of particle density. The dynamics of 
the local temperature is the only quantity we need to model. 

We now return to study the collision terms seen on the RHS of Eq. (|26|). A 
related quantity is the (L-invariant) production rate A 12 ^ 34 of particles per unit 
time and space, defined usually with respect to chemically equilibrated distributions: 

a 12^34 _ I „oo„oo/^ „, \ 12^34 / qq\ 

A = -T- — c Pi P 2 {VsVl2)T ■ {• i2 ) 

1 + 01,2 
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The factor 1/(1 + Si^) is introduced to compensate double-counting of identical 
particle pairs. In terms of the L-invariant A , Eq. ( p6[ ) takes the form: 

+l q (Th q (T)A^ sS ~ 1s (t) 1s (t)(A sS ^9 + a sS ^). (33) 

Only weak interactions convert quark flavors, thus, on hadronic time scale, we 
have 7s !g (r) = 7s.^(t). Moreover, detailed balance, arising from the time reversal 
symmetry of the microscopic reactions, assures that the invariant rates for for- 
ward/backward reactions are the same, specifically 



.4 



12 — 34 



A 



34^12 



(34) 



and thus: 



Tp° s 



dT 



T Z K 2 {z) 



= Y g (r)A^ ss 



1 - 



1 - 



7 s 2 (r) 



7 9 2 (r) 



(35) 



When all ji — * 1, the Boltzmann collision term vanishes, we have reached equilib- 



As discussed, the gluon chemical equilibrium is thought to be reached at high 
temperatures well before the strangeness equilibrates chemically, and thus we as- 
sume this in what follows, and the initial conditions we will study refer to the time 
at which gluons are chemically equilibrated. Setting X g = 1 (and without a signif- 
icant further consequence for what follows, since gluons dominate the production 
rate, also X q = 1), we obtain after a straightforward manipulation the dynamical 
equation describing the evolution of the local phase space occupancy of strangeness: 



2r,T 



_ j ±z K 1 {z)_ 
dT ^ T Z K 2 (z) 



= 1-7* 



(36) 



Here, we defined the relaxation time t s of chemical (strangeness) equilibration as 
the ratio of the equilibrium density that is being approached, with the rate at which 
this occurs: 

_ i pT 



' s 2 (As9^ ss + Aii^ ss + . . .) ' 

The factor 1/2 is introduced by convention in order for the quantity r s to describe 
the exponential approach to equilibrium. 

Eq. ( p6| ) is the final analytical result describing the evolution of phase space 
occupancy. Since one generally expects that "f s — > 1 in a monotonic fashion as 
function of time, it is important to appreciate that this equation allows the range 
7 S > 1: when T drops below m s , and 1/t s becomes small, the dilution term (2nd 
term on LHS) in Eq. ( |36| ) dominates the evolution of j s . In simple terms, the high 
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abundance of strangeness produced at high temperature over-populates the available 
phase space at lower temperature, when the equilibration rate cannot keep up with 
the expansion cooling. This behavior of 7 S has been shown for the SPS conditions 
allowing explosive transverse expansion in subsection 3.3. Since we assume that the 
dynamics of transverse expansion of QGP is similar at RHIC as at SPS, we obtain 
similar behavior for j s in section 5 below. 

r s (T) , Eq. (j37]), has been evaluated using pQCD cross section and employing 
next to leading order running of both the strange quark mass and QCD-coupling 
constant a s EJ\ We believe that this method produces a result for a s that can be 
trusted down to just below 1 GeV energy scale which is here relevant. We employ 
results obtained with a s (M z o) = 0.118 and m s (lGeV) = 200 MeV; we have shown 
results with m s (lGeV) = 220 MeV earlier £3 There is some systematic uncertainty 
due to the appearance of the strange quark mass as a fixed rather than running 
value in both, the chemical equilibrium density in Eq. (|37|), and in the dilution 
term in Eq. (|36|). We use the value m s (lGeV), with the 1 GeV energy scale cho- 
sen to correspond to typical interaction scale in the QGP at temperatures under 
consideration. 



5. Expectations for Strange Hadron Production at RHIC 

We now combine all recent advances in theoretical models of strangeness produc- 
tion and data interpretation at SPS enereies with the objective of making reliable 



predictions for the RHIC energy range.E3 First we address the question how much 
strangeness can be expected at RHIC. The numerical study of Eq. (^) becomes 
possible as soon as we define the temporal evolution of the temperature for RHIC 
conditions. We expect that a global cylindrical expansion should describe the dy- 
namics: aside of the longitudinal flow, we allow the cylinder surface to expand given 
the internal thermal pressure. SPS experience suggests that the transverse matter 
flow will not exceed the sound velocity of relativistic matter v± ~ c/y/3. We re- 
call that for a pure longitudinal expansion local entropy density scales according to 
S cx T 3 oc 1/V.E3 It is likely that the transverse flow of matter will accelerate the 
drop in entropy density. We thus consider the following temporal evolution function 
of the temperature: 



T(t) = T Q 



1 



(1 + r 2c/d)(l + T v ± /R ± ) 2 



1/3 



(38) 



We take the thickness of the initial collision region at To = 0.5 GeV to be d(Tg = 
0.5)/2 = 0.75 fm, and the transverse dimension in nearly central Au-Au collisions 
to be R± = 4.5 fm. The time at which thermal initial conditions are reached is 
assumed to be To = lfm/c. When we vary To, the temperature at which the gluon 
equilibrium is reached, we also scale the longitudinal dimension according to: 

d(T ) = (0.5GeV/T ) 3 1.5fm. (39) 



This assures that when comparing the different evolutions of 7 S we are looking at 
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an initial system that has the same entropy content by adjusting its initial volume 
Vq. The reason we vary the initial temperature To down to 300 MeV, maintaining 
the initial entropy content is to understand how the assumption about the chemical 
equilibrium of gluons, reached by definition at To, impacts strangeness evolution. 
In fact when considering decreasing To (and thus increasing Vo), the thermal pro- 
duction is turned on at a later time in the history of the collision. 

The numerical integration of Eq. (|3^) is started at To , and a range of initial 
temperatures 300 < Tq < 600, varying in steps of 50 MeV. The high limit of 
the temperature we explore exceeds somewhat the 'hot glue scenario', E-3 while the 
lower limitjaf To corresponds to the more conservative estimates of possible initial 
conditions £3 Since the initial p-p collisions also produce strangeness, we take as 
an estimate of initial abundance a common initial value 7 s (To) = 0.15. The time 
evolution in the plasma phase is followed up to the break-up of QGP. This condition 
we establish in view of results of the analysis for SPS presented in section 3. We 
recall that SPS-analysis showed that the system dependent baryon and antibaryon 
m^-slopes of particle spectra are result of differences in collective flow in the de- 
confined QGP source at freeze-out. In consequence there is universality of physical 
properties of hadron chemical freeze-out between different SPS systems. This value 
is nearly applicable to RHIC conditions, as can be seen extrapolating the phase 
boundary curve to the small baryochemical potentials. The QGP break-up temper- 
ature Tj PS ~ 143 ± 5 MeV will see just a minor upward change, and we adopt here 
the value Tf mc ~ 150 ± 5 MeV. 

With the freeze-out condition fixed, one would think that the major remaining 
uncertainty comes from the initial gluon equilibration temperature To, and we now 
study how different values of To influence the final state phase space occupancy. We 
integrate numerically Eq. (f36j) and present 7 S as function of both time t in Fig. ||a, 
and temperature T in Fig.|]b, up to the expected QGP breakup at T^ HIC ~ 150 ±5 
MeV. We see that: 

• widely different initial conditions (with similar initial entropy content) lead to 
rather similar chemical conditions at chemical freeze-out of strangeness, 

• despite a series of conservative assumptions, we find, not only, that strange- 
ness equilibrates, but indeed that the dilution effect allows an overpopulation of the 
strange quark phase space. For a wide range of initial conditions, we obtain a narrow 
band 1.15 > j s (Tf) > 1 . We will in the following, taking into account some contri- 
bution from hadronization of gluons in strange/antistrange quarks, adopt what the 
value j s (T f ) = 1.25. 

We now consider how this relatively large value of j Sl characteristic for the 
underlying QGP formation and evolution of strangeness, impacts the strange baryon 
and anti-baryon observable emerging in hadronization. Remembering that major 
changes compared to SPS should occur in rapidity spectra of mesons, baryons and 
antibaryons, we will apply the same hadronization model that worked in the analysis 
of the SPS data. This hypothesis can be falsified easily, since based and compared 
to the Pb-Pb 158A GeV results it implies: 
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a) shape identity of all RHIC m± and y spectra of antibaryons p , A , S , since 
there is no difference in their production mechanism, and the form of the spectra 
is determined in a similar way to SPS energy range by the local temperature and 
flow velocity vector; 

b) the m^-inverse-slopes of these antibaryons should be very similar to the 
result obtained at CERN for Pb-Pb 158^4 GeV, since the expected 3% increase in 
the freeze-out temperature is accompanied by a comparable increase in collective 
transverse flow. 

The abundances of particles produced from QGP within the sudden freeze-out 
model are controlled by several parameters we addressed earlier: the light quark 
fugacity 1 < X q < 1.1 , value is limited by the expected small ratio between baryons 
and mesons (baryon-poor plasma) when the energy per baryon is above 100 GeV, 
strangeness fugacity X s ~ 1 which value for locally neutral plasma assures that 
(s — s) = 0; the light quark phase space occupancy j q ~ 1.5, overabundance value 
due to gluon fragmentation. Given these narrow ranges of chemical parameters 
and the freeze-out temperature Tf = 150 MeV, we compute the expected particle 
production at break-up. In general, we cannot expect that the absolute numbers 
of particles we find are correct, as we have not modeled the important effect of 
flow in the laboratory frame of reference. However, ratios of hadrons subject to 
similar flow effects (compatible hadrons) can be independent of the detailed final 
state dynamics, as the results seen at SPS suggest, and we will look at such ratios 
more closely. 

Taking 7 9 = 1.5 , we choose the value of A,j, see the header of table |7], for 
which the energy per baryon (E/B) is similar to the collision condition (100 GeV), 
which leads to the range X q — 1.03 ± 0.005. We evaluate for these examples aside of 
E/B, the strangeness per baryon s/B and entropy per baryon S/B as shown in the 
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Table 7. For 7 3 = 1.25, A s = 1 and j q , \ q as shown: Top portion: strangeness per baryon s/B, 
energy per baryon E / B[GeV] and entropy per baryon S/B. Bottom portion: sample of hadron 
ratios expected at RHIC. 





1.25 


1.5 


1.5 


1.5 


1.60 


\ 


1.03 


1.025 


1.03 


1.035 


1.03 


E/B[GeV] 


117 


133 


111 


95 


110 


s/B 


18 


16 


13 


12 


12 


S/B 


630 


698 


583 


501 


571 


p/P 


1.19 


1.15 


1.19 


1.22 


1.19 


A/p 


1.74 


1.47 


1.47 


1.45 


1.35 


A/p 


1.85 


1.54 


1.55 


1.55 


1.44 


A/ A 


0.89 


0.91 


0.89 


0.87 


0.89 


S-/A 


0.19 


0.16 


0.16 


0.16 


0.15 


S-/A 


0.20 


0.17 


0.17 


0.17 


0.16 




0.94 


0.95 


0.94 


0.93 


0.94 


n/s~ 


0.147 


0.123 


0.122 


0.122 


0.115 


n/IF 


0.156 


0.130 


0.130 


0.131 


0.122 


n/n 


1 


1. 


1. 


1. 


1. 




0.15 


0.13 


0.13 


0.13 


0.12 


A+A 


0.19 


0.16 


0.16 


0.16 


0.15 


K+/K- 


1.05 


1.04 


1.05 


1.06 


1.05 



top section of the table [7| We do not enforce (s — s) = exactly, but since baryon 
asymmetry is small, strangeness is balanced to better than 2% in the parameter 
range considered. In the bottom portion of table [?], we present the compatible parti- 
cle abundance ratios, computed according to the procedure developed in section 2. 
We have given, aside of the baryon and antibaryon relative yields, also the relative 
kaon yield, which is also well determined within this approach. 

The meaning of these results can be better appreciated when we assume in an 
example the central rapidity density of direct protons is dp/dy\ ccnt , — 25. In ta- 
ble ||, we present the resulting (anti)baryon abundances. The net baryon density 
db/dy ~ 16 ± 3, there is baryon number transparency. We see that (anti)hyperons 
are indeed more abundant than non-strange (anti)baryons. Taking into account the 
disintegration of strange baryons, we are finding a much greater number of observed 
protons dp/dy\°^ ~ 65 ± 5 in the central rapidity region. It is important when 
quoting results from table || to recall that: 

1) we have chosen arbitrarily the overall normalization in table ^|, only particle 
ratios were computed, and 

2) the rapidity baryon density relation to rapidity proton density is a conse- 
quence of the assumed value of X q , which we chose to get E/B ~ 100 GeV per 
participant. 

The most interesting result seen in table ^|, the hyperon-dominance of the baryon 
yields at RHIC, does not depend on detailed model hypothesis. We have explored 
another set of parameters in our first and preliminary report on this matter ,Ej find- 
ing this result. Another interesting property of the hadronizing hot RHIC matter, 



Observing Quark- Gluon Plasma with Strange Hadrons. 33 



Table 8. dN/dy\ ccn ^ assuming in this example rfp/dy| C ont. = 25 . 
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as seen in table |?], is that strangeness yield per participant is expected to be 13-23 
times greater than seen at present at SPS energies, where we have 0.75 strange 
quark pairs per baryon. As seen in table ^, the baryon rapidity density is in this 
examples similar to the proton rapidity density. 

6. Summary and Conclusions 

We believe that this study of SPS strangeness results decisively shows interesting 
new physics. We see considerable convergence of the results around properties of 
suddenly hadronizing QGPE3 The key results we obtained in the Fermi-2000 model 
data analysis are: 

1) the same hadronization temperature T= 142-144 MeV for very different col- 
lision systems with different hadron spectra; 

2) QGP expected results for the source phase space properties: A s = 1 for both 
S- and Pb-collisions, implying A^ b ~ 1.1 ; 

3) 7^ b > 1, indicating that a high strangeness yield was reached before freeze- 
out; 

4) j q > 1 as would be expected from a high entropy phase and the associated 
value S/B ~ 43 ± 3 ; 

5) the yield of strangeness per baryon s/B ~ 0.7 just as predicted by gluon 
fusion in thermal QGP, a point we studied in detail in section 3.3; 

6) the transverse expansion velocity for Pb-Pb: v^ h < l/y3, just below the 
sound velocity of quark matter. 

The universality of the physical properties at chemical freeze-out for S- and Pb- 
induced reactions points to a common nature of the primordial source of hadronic 
particles. The difference in spectra between the two collision systems considered 
arises, in this analysis, due to the difference in the degree of chemical equilibration 
of light and strange quarks, expected for systems of differing size and lifespan, and 
a difference in the collective surface explosion velocity, u| ~ 0.5 < v^ h ~ l/y/3, 
which for larger system is higher, having more time to develop. Considering how 
small the experimental WA97 spectral slope errors shown in table ^ are presently, 
there is now overwhelming evidence that the production mechanism of both A and 
A is the same, which observation is very probably also true for both 5 and S. This 
symmetry between matter-antimatter production is an important cornerstone of the 
claim that the strange antibaryon data can only be interpreted in terms of direct 
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emission from a deconfincd and thus matter-antimatter symmetric quark matter. 

We note that the QGP break-up temperature we find, Tf — 143 MeV, cor- 
responds to an energy density e = 0(0.5) GeV/fm 3 @ Among other interesting 
results which also verify the consistency of the experimental data understanding 
within the Fermi-2000 model, we recall: 

• the exact balancing of strangeness (s — s) — also in the final hadronic particles 
in the symmetric Pb-Pb case; 

• the increase of the baryochemical potential /i^ b = 203 ± 5 > [Jg = 178 ± 5 MeV 
as the collision system grows; 

• the energy per baryon near to the value expected if energy and baryon number 
deposition in the fireball are similar; 

• hadronization into pions at j q — > 7^ = e m "/ 2T/ ~ 1.6 seen in Pb-Pb reactions, 
which is an effective way to convert excess of entropy in the plasma into hadrons, 
without need for reheating, or a mixed phase; the finding of the maximum allowable 
7 9 is intrinsically consistent with the notion of an explosively disintegrating QGP 
phase. 

A reassuring feature of the Pb-Pb analysis related to chemical equilibration has 
been described in subsection 2.3: we find a pion yield which maximizes the entropy 
density of hadronic particles produced.E3 This detailed technical result explains how 
sudden hadronization can occur: in general the deconfined state with broken color 
bonds and thus the high entropy density has to find an exit into the hadronic world, 
maintaining or increasing the total entropy and preferentially also the local entropy 
contained within a small, comoving volume cell. Our analysis of experimental results 
suggests that this is accomplished by generating an over-saturated pion phase space, 
in which the entropy density rises to values as high as are believed to occur in QGP 
at hadronization. Chemical equilibrium hadronization requires the formation of a 
mixed plasma-hadron gas phase and is generally believed to require a relatively long 
time, followed by kinetic reequilibration. In our opinion such a hadronization model 
is now inconsistent with the experimental strange baryon and antibaryon data on 
yields and spectra. The reader should note that such technical differences between 
different groups about the dynamics of the evolution of the hadron fireball after the 
deconfined phase has hadronized, do not impact the primary agreement about the 
deconfined nature of the high density source of hadronic particles. 

We believe that omission to consider chemical non-equilibrium in the study of 
freeze-out conditions employing the analysis of spectral shape (flow) and also pion 
correlation (HBT) effectisihe source of the difference of results here presented with 
some other recent workEio To understand the source of this difference it is impor- 
tant to realize that there is a considerable influence on the shape of pion spectra 
by the light quark chemical non-equilibrium which the data analysis presented in- 
cludes: the cocktail of resonance decays contributing to pion spectra is altered, and 
moreover, there is spectral deformation at low mj_ due to pion correlation effects 
caused by the overpopulated phase spaceE3 We note that results presented also dif- 
fer somewhat from the WA98 experiment analysis addressing solely ir° spectra, and 
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which again assumes pion chemical equilibriumE3 In consequence, the ^"-freeze- 
out conditions as seen in Ref.El, Tablet arc different from those determined here. 
On the other hand, another recent hadron spectral shape analysis,B3 which did not 
introduce low mi pion spectra into consideration obtains a chemical freeze-out con- 
ditions nearly identical to those we discussed. While the precise understanding of 
hadronization condition is required for a measurement of physical properties of QGP 
including the latent heat, the differences discussed are of little if any consequence 
concerning the fundamental issue, the question if deconfinement is achieved. 

The sudden hadronization of entropy rich QGP leads to value 7 g — > 7°, in order 
to connect the entropy rich deconfined and the confined phases more efficiently. The 
dominant pion contribution to the entropy density (and pressure) is nearly twice as 
high at 7 g ~ 7J than at 7 g = 1. Without this phenomenon one has to introduce 
a mechanism that allows the parameter VT 3 to grow, thus expanding either the 
volume V due to formation of the mixed phase or invoking a rise of T in the re- 
heating. The range of values for 7 g is bounded from above by the Bose distribution 
singularity 7^ — > 7^, but a pion condensate is not formed since it 'consumes' energy 
without consuming entropy of the primordial high entropy QGP phase. An inter- 
esting feature of such a mechanism of phase transition which side-steps the need 
to form a mixed phase or reheating is that the chemical non-equilibrium reduces 
and potentially eliminates any discontinuity in the phase transition. This being 
the case, experimental searches will not find the critical fluctuations expected for a 
discontinuous phase change, even if theory implies a 1st order phase transition for 
the statistical equilibrium system. This is in agreement with the failure of NA49 
experiment to find precritical fluctuations in event-by-event analysis .0 

The only not fully quantitatively described particle yields are Q and ft: for 
the parameters we find, the Fermi-2000 model applied to Pb-Pb reactions under- 
predicts this smallest of all hadronic abundances by about factor 2 . This means 
that we expect, within the sudden hadronization of QGP with uncorrelated strange 
quarks in the deconfined phase, only an enhancement of Q and £1 yields by a factor 
8-10 as compared to what is expected from extrapolation of p-A reactions. Since 
the experiment reports an enhancement by factor 15-20, we need to think again. 
This 'failure' of Fermi-2000 model is in fact confirming the early expectation that 
f2 and f2 yields are the best signature of deconfinement JiJ we just must in future 
address the question what exactly this tells us about QGP structure. We have 
argued nearly 20 years ago that strangeness pairing in thecolor anti-triplet channel 
(ss)g in the QGP source would enhance fl and O yields,M a point that is of some 
topical interest today in context of color superdconductivity studiesJla 

Despite this unexpected additional enhancement, we firmly conclude in view of 
all diverse evidence that (multi)strange hadronic particles seen at CERN-SPS are 
emerging from a deconfined QGP phase of hadronic matter and do not undergo a 
re-equilibration after they have been produced. This finding has encouraged us to 
consider within the same computational scheme the production of strange hadrons 
at RHIC conditions. First, we have shown that one can expect strangeness chem- 
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ical equilibration in nuclear collisions at RHIC if the deconfined QGP is formed. 
There will, as at SPS, be overpopulation effect associated with the early strangen- 
ess abundance freeze-out before hadronization. Most importantly for signatures 
of new physics at RHIC, we found that strange (anti)baryon abundances will be 
greater than the yields of non-strange baryons (protons, neutrons). Consequently, 
the rapidity distributions of (anti)protons are arising from decays of (anti)hyperons. 

We are not aware that microscopic model studies reported in the literature 
about RHIC conditions which have noted this remarkable hyperon dominance re- 
sult, see, e.g., Ref. The reader could wonder why is this unusual phenomenon 
not happening at SPS energies described in section 3? At SPS there is still an 
appreciable relative baryon abundance among all hadrons (about 15%) and the 
strangeness yield is at SPS energies only at a level similar to the baryon yield. 
Thus while abundant (anti)hyperon formation begins to set in, there are still many 
non-strange (anti)baryons produced. With increasing per baryon energy the yield 
of strange quark pairs per baryon rises, and at the same time the relative abun- 
dance of baryons among all hadrons diminishes. As result, at RHIC energies, we 
have predicted that hyperons and/or antihyperons are the dominant population 
fraction among all baryons and/or antibaryons. We thus believe that the prepon- 
derance of hyperons as the dominant (anti)baryon population at RHIC energies can 
be uniquely correlated with the formation and sudden hadronization of deconfined 
QGP phase. 

In a nutshell: we find that strangeness and (anti)baryon QGP signatures are 
conclusively proving formation of deconfined quark matter phase at SPS energies, 
and that these signatures of new physics are much more distinct at higher RHIC 
energies. 
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